Nonlinear effect on the transmission of light in a cavity array 
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Taking nonlinear effect into account, we study theoretically the transmission properties of photons 
in a one-dimensional coupled cavities, the cavity located at the center of the cavity array is coupled 
to a two-level system. By the traditional scattering theory, we calculate the transmission rate of 
photons along the cavities, and discuss the effect of nonlinearity on the photon transport. The 
results show that the controllable two-level system can act as a quantum switch in the coherent 
transport of photons. The dynamics of such a system is also studied by numerical simulations, the 
effect of the atom-field detuning and nonlinearity on the dynamics is shown and discussed. 
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I. INTRODUCTION 

Recent experimental progress in the fabrication of 
microcavity arrays and the realization of the quan- 
tum regime in the coupling of atomic-like structures 
to Quantized electromagnetic modes inside the cavities 
0, 0, S II B H S H open up the possibility of using 
them as quantum simulators of many-body physics. The 
possibility to create such cavity arrays have stimulated 
a number of theoretical investigations on the transport 
physics 0, [ToL [Hj], with particular emphasis on possible 
analogies with mesoscopic phenomena in the electronic 
context. This study is motivated by the preparation of 
new quantum states for photons or combined photonic- 
atomic excitations, which offer advantages over other sys- 
tems for the realization of strongly interacting many- 
body models in quantum optical systems. For Joseph- 
son junction arrays and optical lattices, it is difficult to 
address individual sites due to the small separations be- 
tween neighboring sites, but this is easy for coupled cav- 
ity arrays due to the large separations which are usually 
dozens of micrometers and can therefore be accurately 
accessed by optical frequencies. 

Coupled cavity arrays have several interesting poten- 
tial applications, including quantum information pro- 
cessing and simulations of quantum strongly correlated 
many-body systems. By using a real-space model 
Hamiltonian, it has been shown that the transmis- 
sion of a single photon can be switched on or off in 
one-dimensional waveguide coupled with superconduc- 
tivity quantum bitsQ. Indeed the transmission can be 
switched to any predicted value with a tiny change in the 
magnetic field when the sharper Fano resonance peak is 
employed. For a single photon inside a one-dimensional 
resonator waveguide with a two-level system, a general 
spectral structure in which the reflection and transmis- 
sion beyond the usual Breit-Wigner and Fano line shapes 
is predicted 10]. The atom in a half-waveguide can play 



an intrinsic role of semi-transparent mirror for the single 
photon, leading the system to be flexible under control. 

Nonlinearity may appear in quantum systems due to 
the many-body effects and/or system-environment cou- 
plings. Transport and tunneling properties of nonlin- 
ear system (for example, Bose-Einstein condensates) are 
of considerable current interests, both experimentally 
and theoretically. Quantum tunneling through a bar- 
rier is a paradigm of quantum mechanics and usually 
takes place on a nanoscopic scale, such as in two sup- 
perconductors separated by a thin insulator [l2| and two 
reservoirs of superfluid helium connected by nanoscopic 
apertures (l3l [T3|. Recently, tunneling on a macroscopic 
scale (pm) in two weakly linked Bose-Einstein conden- 
sates in a double- well potential has been observed [l5j. 
Similar to tunneling oscillations in superconducting and 
superfluid Josephson junctions, Josephson oscillations 
are observed when the initial population difference is 
chosen to be below a critical value. When the initial 
population difference exceeds the critical value, an in- 
teresting feature can be observed, i.e., tunneling oscilla- 
tions are suppressed due to the nonlinear condensate self- 
interactions. This phenomenon is known as the macro- 
scopic quantum self-trapping. 

Here we explore a different regime of transport for pho- 
tons in a cavity array coupled to a two-level system. Non- 
linear effects arising from Kerr medium in cavities are 
taken into account. The paper is organized as follows. 
In Sec. II, we present a general formulism for the atom- 
cavity array. By the traditional scattering theory, the 
transmission rate is calculated and discussed in Sec. III. 
In Sec. IV, we present a numerical simulation for the dy- 
namics of the coupled atom-cavity array system. An ex- 
perimentally realizable proposal to observe the prediction 
is suggested in Sec.V. Finally, we conclude our results in 
Sec.VI. 



II. FORMULISM 
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Consider a cavity array coupling to a two-level system. 
Photon hopping can occur between neighboring cavities 
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FIG. 1: (color online) Schematic illustration of the coupled 
cavity arrays, in which the photons propagate. A two-level 
atom is located at the zeroth cavity and coupled to the cavity 
field. 



due to the overlap of the special profile of the cavity 
modes, see figure [TJ Introducing the creation and anni- 
hilation operators of the cavity modes, and a,j, the 
Hamiltonian of the field in the cavity can be written as, 



H l c = uj y~] djdj — ^Sj-(o]Oj+i + h.c). 



(1) 



In Hamiltonian Eq. (fTJ) that all the cavities have the same 
resonant frequency and the same spacing overlap between 
all neighboring cavities have been assumed. In contrast 
to the previous studies, we focus here on the influence of 
nonlinearity on the transparent properties of photons in 
the cavity array. The nonlinearity comes from the optical 
Kerr effects, where the electronic field is due to the light 
itself. Taking this nonlinearity into account, we add an 
additional term g^jia^a) 2 into the Hamiltonian, and 
the Hamiltonian for the cavity array therefore becomes 



H r = Hi 



/ > 
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(2) 



The location of the two-level atom in the cavity array is 
chosen to be the origin of the coordinate axis, which is 
only coupled to the zeroth cavity field in the cavity array. 
We assume that the cavity decay and atomic spontaneous 
emission are ignored. This configuration is different from 
that in RefspJ, [O, d, 13, where several approxima- 
tions have been made to treat the problem. The coupling 
of the two-level atom to zeroth cavity field is described 
by the Janes-Cummings model(under the Rotating Wave 
Approximation) , 



ff J =0|e)(e| + Jo(aS|<7><e| + |e>( 9 |a ), 



(3) 



where Jo denotes the coupling constant. To observe the 
effect of this nonlinearity, we consider the situation where 
the field intensities are very high, this yields 

ihaj — ujaj — £(<Xj+i + ctj-i) + 2g\aj\ 2 aj, 

iha = uja - £(a_i + a x ) + 2g\a Q \ 2 a + J (cr_), 

ih(&-) = n(<r_) - Jqolo{(Tz), 

ih(a z ) = -2J ao(r7_) + 2J a ((j + ), 

ih{&+) = -0(cr + ) + J a (cr 2 ). (4) 

Here, to account for the classical and quantum fluctua- 
tions, each operator was decomposed as a sum of its aver- 
age value and a small fluctuation, i.e., cij — otj +Saj, and 



cr +i _ iZ = ((T +i _ i2 )-|-(5(T+ i _ i2 . Substituting these quantities 
into the Heisenberg equation, ihdA/dt = [A, H c + Hj], 
(A = aj,a+,a-,a z ,j = 0,1,2,...), we have Eq.(j4]) and 
the following equation for the fluctuations, 



dt 3 



at 



ot 

ih—d(j z 
at 

l h-8a + 



w5aj — £(($Oj+i + <5aj_i) 
4g\a J \ 2 6a ] +2ga 2 Sa], {j ^ 0), 
LuSag — £(8ai + 5a,-x) + Jq5<t- 
4g\a \ 2 8a + 2gal5a\, 
QSa- - J a 5<7 z - J 5a (a z ), 

-2Jo(5al(a-) - a Sa + -f 

— f2(5er+ + JochqSctz + Jo8a\{a 



Sa (a + )), 
(5) 



In deriving Eq.© we have eliminated the average value 
contribution and linearized the fluctuations. The dy- 
namical stability of the system can be obtained by 
Eq.®, this method adopts the linear stability anal- 
ysis that has wide applications in various nonlinear 
systems. Eq.© can be written as ih^r — H e ffV 
with V = (Sa-N , 5aJ_ N , Sa-N+i, Sa_ N+v 5a , 5%, 
Sajf-i, Sa N _ lt SdN, So, N ,S(T^,Sa z , S(T + ) T , and-ff e // can 
be derived from Eq.©. We shall not present the details 
of analysis here, which will be discussed elsewhere. 



III. TRANSMISSION RATE 

Now we turn to the shifted picture defined by the trans- 
formation, (er_) — > (er~_) = ((7_)e J °' t *, fi/. = — 2£cosfc + 
oj + 2g\aj\ 2 , by setting d{<f-)/dt = 0, and solve for (crl) 
in terms of (<7 Z ), we arrive at 

n k a = ua Q - £(a_i + a{) + 2g\a Q \ 2 a a + a J /(A, J ), 



where /(A, J ) = q°1q^ , and A = Q — fl)-. Equations 
([6]) represent a coupled set of stationary states of the 
total system for the fields in the high-intensity limit. The 
scattering equation for j ^ has the solution, 
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j<0 
3>0, 



(7) 



where s and r denote the transmission and reflection am- 
plitude, respectively. By solving the scattering equation 
with the continuous condition o^(fc) = ag(fc), we obtain 
the transmission amplitude s satisfying, 



2g\s\ 2 ■ s + (2£cosfc 
+2i£sinfe = 0, 



J 



Q — Qi 



- 2£e lte )s 



(8) 
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FIG. 2: (color online) The transmission coefficient \s\ 2 as a 
function of the momentum k of the incident photons, (a)- 
(d) correspond to different nonlinear coupling constant g = 
0, O.f, —0.5, 0.5, —1, and 1, respectively. The other parameters 
chosen are J = 1,£ = 1, u = 2, and = 3. 



where J = Jq y/ (d z ). Set s — x + iy with real numbers x 
and y, Eq.© follows, 

2.9(£ 2 + y 2 )y + ( Q _ n )y - 2x£ sin fc + 2£ sin fc = 0, 

2 5 (a; 2 +y 2 ) a : + (- -)x + 2y£sinfc = 0, (9) 

life — i i 

Eq.© are coupled cubic equations giving three roots to 
the transmission amplitude s. Considering the restric- 
tions on the amplitude s (i.e., x,y are real number and 
x 2 + y 2 < 1), we have performed extensive numerical 
calculations for the transmission amplitude s, selected 
results are presented in figure 2-5. 

Figure [2] shows the transmission coefficient \s\ 2 as 
a function of the momentum k of the incident pho- 
tons. Figure [2J-(a) illustrates the transmission coeffi- 
cient \s\ 2 without nonlinear interactions, while figure [2} 
(b),(c) and (d) show the effect of nonlinear coupling on 
the transmission coefficient |s| 2 with coupling constant 
g = 0.1,-0.5,0.5,-1 and 1, respectively. We can find 
from [2]-(b),(c) and (d) that the lines are not continu- 
ous, indicating no solutions to Eq.© can be found at 
the blank points. This is different from the case without 
nonlinear couplings (namely, g = 0). Figure [5] also shows 
that the transmission coefficients not only depend on the 
strength of the nonlinear coupling, but also on the fea- 
ture of the nonlinear interaction, i.e., the transmission 
coefficients behave different for attractive and repulsive 
nonlinear interactions (see, (c) and (d), (e) and (f)). 

For large nonlinear coupling constants, the transmis- 
sion coefficient |s| 2 is zero, as figure [3] shows. This means 
that the photons can not transmitted along the cavity ar- 
ray. For g — and J = 0, the transmission coefficient |s| 2 
is 1, however this is not the case for nonzero g and J = 0, 
|s| 2 would depend on k. For small g, the nonlinear effects 
block the photon transport in most cases, however, for 
some special fc, the nonlinearity favors the photon trans- 




FIG. 3: (color online) The transmission coefficient |s| 2 as a 
function of the momentum k and the nonlinear coupling con- 
stant g. J = 1, £ = 1, u) = 2, and = 3 were chosen for this 
plot. 
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FIG. 4: (color online)This contour plot shows the transmis- 
sion coefficient as a function of k and g with different cavity- 
atom coupling constant J. (a) J — 0, (b) J = 0.5, (c) J = 1, 
(d) J — 1.5, (e) J — 2, and (f) J = 5. The other parameters 
chosen are £ = 1, w = 2, and = 3. Note that the blank 
regions in the figure indicate no solutions have found for |s| . 



mission. For large <?, the nonlinear effect always decreases 
the transmission rate independent of what fc takes. The 
coupling of the atom to the photons can shift the peaks 
of the transmission spectrum, this is shown in figure |U 
For fixed fc and g 7 the dependence of the transmission 
coefficient |s| 2 on J and £ is shown in figure Figured] 
shows that the peak of the transmission spectrum is sym- 
metric about g — when J — 0. As the coupling constat 
J increases, the peaks are shifted and lose the symmetry 
about g, i.e., the stronger the atom-field coupling, the 
bigger the peak shift. Two observations can be found 
from figure (1) For a fixed J, the transmission coeffi- 
cient |s| 2 increases with £. This can be easily understood 
as that large £ would increase the photon transmitting, 
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FIG. 5: (color online) Transmission coefficient as a function of 
J and £ with fc = 1, g = l,w = 2, and f2 = 3. 
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resulting in larger |s| 2 . (2) For a fixed £, the dependence 
of the transmission coefficient |s| 2 on J behaves different 
for positive and negative £. For £ > 0, |s| 2 increases as 
| J\ decreases, whereas for £ < 0, |s| 2 becomes larger with 
| J | decreases. 



IV. DYNAMICS 

To gain insight of the nonlinear effect on the transport 
of photons along the cavity array, we numerically sim- 
ulate the dynamics of the system, and present the nu- 
merical results in this section. For this purpose, we first 
rescall the Hamiltonian H — H c +H] by the total number 
of photons in these cavities. Notice that for a very large 
number of photons, M — (M) ^> 1 with M — J2i a l a ii 
the total number of photons is (approximately) a con- 
stant of motion for the system. So, we may simulate the 
dynamics by using a Hamiltonian obtained by replacing 
„t 



|,Oj with al/vM,ai/vM, respectively, this leads to a 
set of coupled equations which have the same form as 
in Eq.((4]), but with rescaled constants, uj r — uj, £ r = £, 
g r = M-g, fi r = fi, and J r = JJyM. We plot the photon 
number in the —1th (blue and dashed line) and +lth (red 
and solid line) cavities as a function of time in Fig[6]and 
FigEl The initial state of photons is all photons (total 
number M = 15) in the —1th cavity, while the atom ini- 
tially is in its excited state \e) . The —1th cavity is the first 
cavity left of the atom, and the +lth one is the first on the 
right, the atom occupies the 0th cavity, as we described 
before. Fig[5] shows that for resonantly atom-field cou- 
pling (Fig[6]-(a)), the photon in the both cavities arrive in 
a synchronizely changed state for a very short time. As 
the detuning (u> — Q) increases, it needs a longer time to 
evolve into such a state (see Fig|B]-(b),(c)). We also find 
from Fig[6] that for resonantly atom-field coupling, it is 
easier for the photons to transfer through the 0th cav- 
ity than the case of non-resonant coupling. This can be 
found by observing the first peak of the red-solid curve, 
obviously, the first peak of the red-solid line in Fig[6]-(a) 
is higher than that in (b) and (c). Fig|7] shows the non- 
linear effect on the dynamics of the system. Again, the 
nonlinearity affects the time needed for the cavity field to 



FIG. 6: (color online) Photon number in —1th cavity (blue 
and dashed line) and -l-lth cavity (red and solid line) as a 
function of time. The parameters chosen are uo r = 2, J r — 15, 
M = 15, £ r = 1, g r = 2. (a) Sl r = 2, (b)O r = 3, and (c) 

a- = 5. 




FIG. 7: (color online) The same as Fig|6j but with different 
parameters. ui r — 2, S!,- = 3. (a) g r = 0.5, (b) g r — 2, (c) 
g r = 2.9, (d) g r = 3. (e) is a long time plot for (d). 



evolve into a synchronized state. The larger the nonlin- 
ear coupling constant g is, the longer the time to arrive at 
such a state. Fig[7| also shows that the nonlinear effect 
blocks the photon transfer along the cavity array with 
strong nonlinearity, reminiscent of the self-trapping for 
Bose-Eintein condensates in a double-well potential. 



V. POSSIBLE EXPERIMENTAL OBSERVATION 

In this section, we propose an experimentally acces- 
sible quantum system to demonstrate our theoretical 
predictions for the transmission coefficient |s| 2 . Con- 
sider a system of (27V + 1) single- mode coupled non- 
linear waveguides ( denoted as —N,—N + 1,—N + 
2, —1, 0, 1, N — 1, N, see figure|S]). The waveguides 
are identical in shape and have a constant width[20] along 
the propagation direction z. The distances between the 
waveguides equal to each other, as a consequence, the 
coupling constants between the waveguides remain un- 
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FIG. 8: (color online) A schematic illustration of the setup. 
The relative distance between the coupled waveguides equal 
to each other, resulting in the same coupling constant among 
the neighboring waveguides. The two-level atom was put in 
the zeroth waveguide in the center. 

changed along the propagation. L denotes the waveg- 
uide length. We assume that only couplings between 
neighboring waveguides are practically nonzero in this 
configuration. The two-level atom is put in the zeroth 
waveguide. The waveguides can be designed to have a 
width of 3/im and length of 18mm. The edge-to-edge 
distance between waveguides is 5/zm, yielding a coupling 
constant of 2500 m _1 . The evolution of the modal am- 
plitudes in those waveguides can be described by Eq.Q 
in the intensive light limit. The nonlinear couplings g 
stem from the Kerr nonlinear coefficient of the waveg- 
uides. These couplings are important only in the nonlin- 
ear regime and can be neglected at low light power levels. 
We would like to note that A in Eq.([5]) should be under- 
stood as dA/dz (A = a,j, o + , er_, <r z ,j = 0, 1, 2, ...) in this 
coupled waveguide system. The waveguides used in this 
proposal can be fabricated on an AlGaAs substrate, us- 
ing standard photolithography techniques[2l|. Indeed by 
this configuration in three coupled waveguides, the ef- 
fect of nonlinearity on adiabatic evolution of of light was 



observed [2Cj. Besides, these features that result from the 
effect of nonlinearity can also hopefully be observed in 
coupled superconducting transmission line resonators. 



VI. CONCLUSION 

To sum up, taking nonlinear couplings of photons into 
account, we have studied theoretically the transmission 
properties of intensive light in a one-dimensional waveg- 
uide. The transmission coefficient was calculated and 
discussed. These results show that the nonlinear cou- 
plings sharply influence the transmission, and the atom- 
field couplings shift the peaks in the transmission spec- 
trum. The dependence of the transmission coefficient on 
the couplings between the neighboring sites is also cal- 
culated and analyzed. Besides, we numerically simulate 
the dynamics of the system. The following features are 
found. (1) Resonant atom-field coupling favors the pho- 
ton transport, the larger the detuning is, the longer the 
time needed to arrive at a synchronized state. (2) Non- 
linearity blocks the transfer of photons among the cavity 
array, and postpone the time to arrive at a synchronized 
state. (3) Both a off-resonant atom-field coupling and a 
strong nonlinearity make the photon transport along the 
cavity array difficult with respect to the case of resonant 
atom-field coupling and without (or small) nonlinearity. 

This work was supported by NSF of China under 
Grant No. 10775023. 
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